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We present rigorous results for quantum systems with both bound and continuum states subjected 
to an arbitrary strength time-periodic field. We prove that the wave function takes the form ip(x, t) = 

Y)j 2~2k=o 1 t ke ~ lrTitc 3 4>j,k(3:,t) + ipd{x,t), with 4>j,k{x,t) a set of time-periodic resonant states with 
quasi-energies Oj — Ej — iTj/2, with Ej the Stark-shifted energy and Fj the ionization rate, and Nj 
the multiplicity of each resonance. tpd{x,t) is the dispersive part of the the solution, and is given 
by a power series in t~ x ^ . Generally, Fj > for each j leading to ionization of the atom, but we 
also give examples where Fj — and implying the existence of a time-periodic Floquet bound state. 
The quantity Oj — Ej — iTj/2 has a convergent perturbation expansion for small field strengths. 

PACS numbers: 02.30. Jr,03.65.-w,32.80.Fb,32.80.Rm 



The ionization of atoms subjected to external time- 
dependent fields is an issue of central importance in 
atomic physics. There exist a variety of methods for 
treating this problem, including perturbation theory 
(Fermi's golden rule), numerical integration of the time- 
dependent Schrodinger equation, semi-classical phase- 
space analysis, Floquet theory and Complex dilations 
[J, H, H, 0, IB IE 0] • Still, there are very few rigorous re- 
sults proving or disproving ionization by a periodic field 
of arbitrary strength and frequency, even for the sim- 
plest systems with both bound and continuum states. 
Such results are clearly desirable from both a theoretical 
and practical point of view. Numerical results are also 
difficult since derealization of the wavefunction of the 
electron creates truncation errors at the boundary of the 
computational domain (absorbing potentials can help, 
but may also cause difficult to detect errors d, Section 
3.2]). This has motivated us to undertake a systematic 
study of this problem in the context of non-relativistic 
quantum mechanics where the radiation field causing the 
ionization is treated classically Q- We present here rig- 
orous results applicable to the behavior of such systems 
in very intense laser fields, a subject of much current in- 
terest 0,|1]. 

For weak fields our results agree with those of pertur- 
bation theory. They show, in fact, that certain quantities 
such as the ionization rate and Stark shift have conver- 
gent expansions in the field strength, thus rigorously jus- 
tifying the use of perturbation theory for weak fields. For 
larger fields, however, the behavior is very complex even 
in the simplest model systems and the results can be 
qualitatively different from those given by perturbation 
theory. In particular the ionization probability may not 
be monotone in the strength of the field, leading to some 
kind of stabilization @, M, EL H • 

In our work, we study the long time behavior of the 
solution of the Schrodinger equation in d dimensions 



--A + V (x) + V 1 (x,t) 



The units are chosen so that H = m = 1. Here, i e 
t > 0, Vo(x) is a binding potential having bound and con- 
tinuum states, and V\{x,t) = YlJLi [Qjixje 1 ^ 1 + c.c] is 
a time-periodic potential representing the radiation field. 
Some of our results require that Qj (x) = for j > M > 1. 
The initial condition t/j(x,0) = i/jq(x) is taken to be lo- 
calized and of unit norm; for convenience we also assume 
ipo (x) = for | a; | > R with R arbitrary. 

Of primary interest is whether the system ionizes un- 
der the influence of the forcing V\{x,t), as well as the 
rate of ionization if it occurs. Ionization corresponds to 
derealization of the wavefunction as t — > oo. In particu- 
lar we say that the system, e.g. an atom, will fully ionize 
if the probability of finding the particle in any bounded 
spatial region B C M. d goes to zero as time becomes large, 
i.e. J B \ip(x,t)\ 2 dx — > as t — > oo. Our results however 
go beyond this and give in many cases a comprehensive 
qualitative picture of the time evolution of ip(x, t). 

The time evolution is described by a resonance expan- 
sions for tp[x,t). That is, under the simplifying assump- 
tion that tp(x,0) is compactly supported, we obtain the 
decomposition of i/j(x, t) 



M Nj-l 
j=0 k=0 
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where M is the number of bound states and resonances 
(possibly infinite), and Nj is the multiplicity of the res- 
onance <Tj. The exponents <jj = Ej — iFj/2 are complex 
numbers, with Ej = -Rcr, interpreted as a Stark shifted 



energy and T^ 



-2Ser 7 - the ionization rate. The func- 



tions 4>j t k{x, i) are time-periodic Gamow vectors [l4| with 
period 2ir/uj. This means the <fij^k(x,t) are smooth eigen- 
values of the Floquet Hamiltonian, 



-id t --A + V {x) + V 1 (x,t) 



4>j,k = 4>j,k- 



yb{x,t) (1) 



(3) 

where 4>j,-i = 0. Eq. |3]) is supplemented with an 
"outgoing wave" boundary condition which depends on 
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the potentials Vq(x) and V%(x, t). When S&Xj < 0, this 
boundary condition causes 4>j t k(x,t) to be exponentially 
growing in space. On the other hand, if SJcrj = 0, then 
4>j,k(x,t) = for k > 1 and 4>j_o(x,t) is a standard 
L 2 eigenvalue of the Floquct Hamiltonian. Our analy- 
sis can be used to make rigorous calculations involving 
these Gamow vectors, e.g. [I3.[TH|, 

In the case when Vq(x) + Vi(x,t) = for \x\ > R in 
three dimensions, then 4>j^(x,t) takes the form 



j',oOM) = E/ 
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for |x| > R. The square root is chosen to make the 

smaller angle with the real axis, and '(z) is the spher- 
ical Bcssel function of the third kind. 

The last term, ipd(x,t) is the dispersive part of the 
solution and has the expansion 



ip d (x,t) ~ ^d k (x,t)r 



■k/2 



(5) 
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Generically d\{x,t) = 0, but the presence of a zero en- 
ergy resonance causes di(x,t) to be nonzero. In the case 
of finite range V${x) and Vi(x, i) = 0, we can resum this 
series using the Borel summability procedure [If 
and we conjecture that this can be done in general 
It follows from (J5J) that for very large times the decay 
of any initially localized states will be given by a power 
law, which is expected since solutions to the Schrodinger 
equation with Vq(x) = Vy{x,t) — behave this way as 
well. In the cases studied the power law decay is domi- 
nant only for times when the survival probability in any 
bound state is already very small, but it has in fact been 
observed experimentally [191 ] . 

Equation (|2|) provides a rigorous definition of the ion- 
ization rate and Stark shifted energy regardless of the 
strength of the field. Even in the case when Vi(x, t) = 0, 
the Hamiltonian Ho = — (1/2) A + Vo(a;) will have many 
resonances [20] (complex energies Ej — iTj/2 for which 
Gamow vectors can be found). As Vi(x, t) is "switched 
on" , the bound states will (generically) become reso- 
nances, while the already existing resonances will change 
as well. In general, letting the parameter e represent the 
strength of the field, the quasi-energies <Jj(e) will vary an- 
alytically in e, except for an isolated set of points where 
(7j(e) may have singularities. The bound states and reso- 
nances change analytically as well, and bound states be- 
come resonances when Tj(e) = 2Q<jj(e) becomes nonzero. 

When the initial state is not compactly supported 
there may be additional exponentially or polynomi- 
ally decaying terms which depend explicitly on ipo(x)- 
These terms are not resonances, and are present even in 
the solution of the free Schrodinger equation (Vq(x) — 
Vi(x,t) = 0), and can be distinguished from resonances 
by this fact. 



Given the expansion ([2]), it follows that the only way 
ionization can fail is if there exists an <jj,<j)jfi(x,t) pair 
solving with 3<jj = 0. Such a pair represents an eigen- 
value/eigenvector pair associated to the Floquet Equa- 
tion ([3]) (with Nj = 1) with cfij^x, t) being 27r/w-periodic 
in t and square- integrable in space. If <pjfi{x,i) were not 
square-integrable in space, then conservation of probabil- 
ity would be violated. Proving ionization then involves 
showing that nonzero solutions of ([3]) do not exist. 

Actually this dichotomy, cither the system ionizes or 
there exists a square-integrable Floquet eigenstate, holds 
in more generality than @. It essentially rules out the 
existence of a singular spectrum of the Floquet Hamil- 
tonian or solutions which are localized but not 2it/uj- 
periodic. This is proven to hold in all the cases we 
have investigated and we expect it to hold universally, 
c.f. [1, 0|- We further expect that the first alterna- 
tive, or Tj > 0, will hold for essentially all Vo(x,t) and 
V\{x,t). As shown below however, there are explicit ex- 
amples where there are square-integrable solutions of 
and thus time-periodic bound states cf)jfl(x,t) (of period 
2tt/lo). In these cases the system never ionizes, unless of 
course ip(x, 0) is orthogonal to <f>j,o{ x , 0). 

We describe more explicitly our results for several sys- 
tems. 

• We considered the widely studied [J [H, [2l[ one di- 
mensional system with Vq(x) = — 25(x) and dipolc 
coupling to an oscillating electric field: 



Vi(x,t) = £{t)x 



C.C. X 



(6) 



In [Tl[ we prove for all field strength and all u> that 
© holds. We also show that if £{t) is a trigono- 
metric polynomial, i.e. £j — for j > M > 1, 
then ionization occurs. For more general £(t) all 
we know is that |2]) holds implying that either the 
system ionizes or else it has a Floquet bound state. 

Taking the same system with Vq(x) — — 25(x), 
but with Vi(x,t) = S(x)£(t), we showed that the 
system ionizes if £{t) is a trigonometric polyno- 
mial. In this case, we also proved that if £(t) = 
2a\(X — cos(wt))/(l + A 2 — 2Acos(w£)), then there 
is a continuous set of values in the (a, A, oj) param- 
eter space for which ionization does not hold and a 
Floquet bound state exists, a strictly nonperturba- 
tive situation [22]. A similar situation occurs [Tol ] 
when Vi(x,t) = £ sm(ujt)[S(x — a) — 5(x + a)] for 
a continuous set of (£, a, u>) even when Vg(x) = 0. 
An example of a continuous Vq(x) and V\(x,t) for 
which there are non-ionizing Floquet bound states 
can be constructed from breather modes of the non- 
linear Schrodinger equation [23j. 

These examples indicate that although uncommon, 
stabilization can occur when a resonance reaches 
the real line. It should be noted that in all these 
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cases, the Floquet bound states are time-dependent 
and have many Fourier components in time. Even 
in the absence of bound states, it is possible for res- 
onances (Tj to move very close to the real line and 
survive for a long time. In fact we believe that 
such occurences correspond to what are usually 
termed Laser Induced Continuum States (LICS). 
LICS is a phenomenon in which, under the influ- 
ence of a strong laser field, an atom with ground 
state Eq < appears to form a bound state at en- 
ergy E + E s + lu > (with E s the Stark shift), 
which is located in the continuum [24l l25l. l26l |27| . 
Our mode of analysis may form a means of studying 
this phenomenon. 

• Our last example is that of a Coulomb binding po- 
tential in 3 dimensions: Vq(x) = — Ze 2 /\x\, and 
Vi(x,t) periodic in time and compactly supported 
in space (i.e. Vi(x,t) = for |x| > R). We show 
in this case that the dichotomy holds as well, i.e. 
either Floquet bound statesexist or ionization oc- 
curs. We prove the absence of Floquet bound states 
when Vx(x,t) — Q(\x\) cos(tut) with ri(|x|)> for 
|x| < R, and n(\x\) = for |x| > R 9}. This 
case is mathematically more difficult than the other 
cases due to the long range of the Coulomb poten- 
tial which causes an accumulation of eigenvalues or 
resonances near zero energies. In this case, we have 
not yet proven ([2]) , though we are optimistic about 
proving a modified version of it (Eq. ^ must al- 
tered to account for the long range potential). 

Main steps of the proofs: The basic idea of our ap- 
proach (which must be adapted to each case) is to relate 
the time asymptotics of %j){x,t) in the time domain with 
the analyticity properties of tp(x, a) in the energy do- 
main. That is, we let ip{x,p) = J Q e~ pt il;(x,t)dt be the 

Laplace transform of ip(x,t), and consider ip(x,a) with 
(7 = —ip so that cr can be interpreted as a quasi-energy. 
The Schrodinger equation then becomes an inhomoge- 
neous integral equation for tj)(x,a). We then show that 
this equation is of a modified Fredholm type in a suitable 
Hilbert space. Application of the Fredholm Alternative 
yields the above mentioned dichotomy: either the homo- 
geneous equation ([3]) has a nonzero solution (leading to 
a Floquet bound state) or else ip{ x : cr) is analytic in the 
a variable on the inverse Laplace contour, except for a 
branch point near a = 0. The Riemann-Lebesgue lemma 
(stating that smoothness in a implies decay in t) implies 
ionization. 

More quantitative results such as ^ are obtained by 
pushing the inverse Fourier contour into the lower half 
plane in a. Resonances appear naturally as poles of 
ip{x,o) with So- < and bound states occur as poles 
with 3cr = 0. Poles occur only where ([3]) has a solution 
satisfying an outgoing wave condition. Standard Fred- 
holm theory shows solutions are analytic in both cr and 
the field strength. The implicit function theorem there- 



fore shows that the pole position is analytic (except for 
isolated singularities) in the field strength, which justifies 
perturbation theory. 

In the case of the Coulomb potential, the proof of com- 
pactness and analyticity is made difficult by the slow de- 
cay of Coulomb potential at oo which affects the decay 
rate of ip{ x : cr) as cr ^ ±oo on or close to the imaginary 
axis. Since the discrete spectrum has an accumulation 
point at zero, the Fredholm operator still has an essential 
singularity at cr = since V\(x, t) affects distant orbitals 
only minimally. Nonetheless, by appropriate choice of 
variables, it is possible to show boundedness as a = 
is approached from the upper half-plane. The bounded- 
ness and continuity are reflected in the form of tp(x, a) 
as well: if indeed the associated homogeneous problem 
has only the zero solution, then ionization follows from 
Riemann-Lebesgue Lemma. 

Perturbation theory in this formulation was carried out 
for the systems with a <5-function binding potential and 
Vi(x,t) — £(t)5(x) or V\{x,t) = £{t) ■ x. In particular 
we showed that the leading order behavior is V ~ (£ 2 ) N , 
where £ 2 is the intensity of the radiation field and N is 
the smallest integer for which Nut > — Eq with Eq = — 1 
the ground state energy of the unperturbed system. In 
the event that Nui — Eq for some N, then the pertur- 
bation expansion is done in powers of £ rather than £ 2 . 
The integer N has the interpretation as the "number of 
photons" required for ionization. In the case when co > I 
(corresponding to N = 1) the lowest order term in the 
expansion reduces to Fermi's Golden Rule [2|. We expect 
this to be true generally. 

To prove that ionization occurs, we need only prove 
that there are no localized solutions of ([31) ■ The method 
we have used in [Tl[ is to solve ([3]) for large x and show 
that if the solution decays near x = oo, then the solution 
must be discontinuous somewhere (e.g. at x = 0). 

An example illustrating this method is given by 
Vq(x) — —25(x) and V\{x,t) = £{t)x. In this case, we 
can avoid the use of WKB for large x in lieu of an exact 
solution, and we use patching at x = instead of asymp- 
totic matching. Elementary calculations show that if we 
ignore Vo(x) = — 26(x), has the family of solutions: 

X m ,±(z,t) = e ±^ e -^t e ±x mC (t) (7a) 
A m = —i\Jcrj + rnuj (7b) 

^— \{inujy J 

If ([3]) has a solution, to the left and right of x = 0, we 
can expand it in terms of Xm,±(x,t): 

M^) = \$ m<M °i% Xm ' + t> t } H° n (8) 

The coefficients <p^ R must be such that (j>o(x,t) is con- 
tinuous at x — and the first derivative in x must have 
an appropriate jump at x = 0. In this expression, Mq 
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is the smallest integer so that A m is real and positive 
for all m < M . The terms </>^ with m > M must 
be zero because any term with m > Mq would oscil- 
late rather than decay for large negative x, thus making 
<po(x,t) non-integrable in space. Any term of the form 
Xm,-(x,t) would grow exponentially as x — > — oo, simi- 
larly making (f>o(x, t) non-integrable. The same argument 
applies to x > 0. 

Physically, this means the following. If the resonant 
state 4>jfl(x,t) included terms with m > M , it would 
be oscillatory near x — oo. This would correspond to a 
state which radiates mass away without decreasing. This 
would violate the law of the conservation of probability, 
and is therefore impossible. 

To prove ionization, we show that no function with the 
expansion ([H]) can be continuous at x = 0. In we use 
a Phragmen-Lindelof argument on the analytic continu- 
ation (in t) of 0o (0, t). The argument shows that (using 
© for x < 0) (f>Q (0,t) must be analytic and bounded for 
< 3ft < tt, St > and tt < 3ft < 2vr, St < 0. Similarly, 
we use the expansion for x > to show boundedness for 
< 3ft < tt, St < and tt < 3ft < 2vr, St > 0. The 
Phragmen-Lindelof theorem can then be used to show 
that 4>o(0, i) is bounded when 3ft = 0, tt, 2tt. This implies 
that <po(0,t) is analytic and bounded, an impossible sit- 
uation unless <po(0,t) = 0. This means that there exist 
no solutions to thereby implying that all resonances 



in ([2]) are decaying. 

For a Coulombic binding potential in 3 dimensions, 
V (x) = —Z/\x\, and Vi(x,t) = f2(|x|) cos(wt) (with 
n(\x\) = for |x| > R and Q(\x\) > for \x\ < R 
the analysis is accomplished similarly. Instead of trying 
to match solutions (J7J) at x = as in (JSJ), we attempt 
to match Coulomb wave functions [23, Chapter 14] for 
x > R to time-periodic solutions of the full Schrodinger 
equation on the region |x < R. The proof of ionization 
again involves showing that no solution can be continuous 
and differentiable at the origin. Again, a crucial ingre- 
dient in the proof is the fact that the positive Fourier 
coefficients of (f)j$(x,t) must be zero for |x| > R, and 
required developing new asymptotic tools for the WKB 
analysis of infinite differential-difference systems. 



Acknowledgments 

We thank R. Costin and A. Rokhlenko who partici- 
pated in some of this work. We also thank K. Yajima, 
A. Soffer and M. Kiessling for helpful comments. Re- 
search supported by NSF Grants DMS-0100495, DMS- 
0406193, DMS-0600369, DMS01-00490, DMR 01-279-26 
and AFOSR grant AF-FA9550-04. 



[1] S. C. R.M. Susskind and E. Valeo, Phys. Rev. A 42, 3090 
(1990). 

[2] J. D.-R. C. Cohen- Tannoudji and G. Gryndberg, Atom- 
Photon Interactions (Wiley - New York, 1992). 

[3] E. A. S. Darko Dimitrovski and J. S. Briggs, Phys. Rev. 
Let. 93 (2004), ISSN 0031-9007. 

[4] E. A. S. Darko Dimitrovski, T. P. Grozdanov and J. S. 
Briggs, J. Phys. B. 36, 1351 (2003), ISSN 0953-4075. 

[5] S. Gram and K. Yajima, Commun. Math. Phys. 89, 277 
(1983). 

[6] A. J. A. Galtbayar and K. Yajima, J. Stat. Phys. 116, 
231 (2004). 

[7] K. Yajima, Commun. Math. Phys. 87, 331 (1982). 

[8] A. Soffer and C. Stucchio, Journal of Computational 

Physics (2007), math.NA/0609183. 
[9] O. Costin, J. Lebowitz, and S. Tanveer (2006), 
math.AP/0611818. 
[10] A. Rokhlenko, O. Costin, and J. L. Lebowitz, J. Phys. A 

35, 8943 (2002), ISSN 0305-4470. 
[11] C. Stucchio, O. Costin, and J. L. Lebowitz (????), sub- 
mitted. 

[12] A. Fring, V. Kostrykin, and R. Schrader, J. Phys. A 30, 

8599 (1997), ISSN 0305-4470. 
[13] A. F. C. Figueira de Morisson Faria and R. Schrader, 

arxiv:physics/9808047 (1998). 
[14] G. Gamow, Z. Phys 51, 204 (1928). 

[15] C.-L. Ho and C.-C. Lee, Physical Review A (Atomic, 
Molecular, and Optical Physics) 71, 012102 (pages 6) 
(2005). 

[16] O. Costin, Duke Math. J. 93, 289 (1998), ISSN 0012- 



7094. 

[17] O. Costin and R. D. Costin, Invent. Math. 145, 425 

(2001), ISSN 0020-9910. 
[18] O. Costin, R. D. Costin, and J. L. Lebowitz, J. Statist. 

Phys. 116, 283 (2004), ISSN 0022-4715. 
[19] C. Rothe, S. I. Hintschich, and A. P. Monkman, Physical 

Review Letters 96, 163601 (pages 4) (2006). 
[20] M. Zworski, Notices Amer. Math. Soc. 46, 319 (1999), 

ISSN 0002-9920. 
[21] T. Mercouris and C. A. Nicolaides, J. Phys. B 32, 2371 

(1999). 

[22] O. Costin, R. D. Costin, J. L. Lebowitz, and 
A. Rokhlenko, Comm. Math. Phys. 221, 1 (2001), ISSN 
0010-3616. 

[23] P. D. Miller, A. Soffer, and M. I. Weinstein, Nonlinearity 

13, 507 (2000), ISSN 0951-7715. 
[24] L. Y. T. Peters and T. Halfmann, Phys. Rev. Let. 95 

(2005), ISSN 0031-9007. 
[25] K. Bohmer, T. Halfmann, L. Yatsenko, D. Charalam- 

bidis, A. Horsmans, and K. Bergmann, Phys. Rev. A. 

66, 1050 (2002). 
[26] T. Halfmann, L. P. Yatsenko, M. Shapiro, B. W. Shore, 

and K. Bergmann, Phys. Rev. A. 58, R46 (1998). 
[27] K. B. T. H. LP. Yatsenko, R.G. Unanyan and B. Shore, 

Opt. Comm. 135, 405 (1997). 
[28] M. Abramawitz and I. Stegun, Handbook of Mathematical 

Functions (Dover, 1965). 
[29] M. Correggi, G. Dell' Antonio, R. Figari, and A. Mantile, 

Comm. Math. Phys. 257, 169 (2005), ISSN 0010-3616. 
[30] M. Lewenstein, P. Balcou, M. Y. Ivanov, A. L'Huillier, 



■5 



and P. B. Corkum, Phys. Rev. A 49, 2117 (1994). 
[31] S. Geltman, J. Phys. B: At. Mol. Phys. 10, 831 (1974). 



